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ABSTRACT 


VI 11 


Seismic response analyses in case of scenario earthquake approach or for 
balanced risk can be conveniently earned out if power spectral density function 
(PSDF) for ground acceleration is estimated with a given level of confidence in 
terms of known parameters like earthquake magmtude, epicentral distance, focal 
depth, site conditions, etc As a first step in tins direction and in extension of 
available models for Fourier spectrum, this study proposes two empirical models 
for scaling PSDF of horizontal ground acceleration The first model considers 
earthquake magmtude, epicentral distance, focal depth, and strong motion du- 
ration as the governing parameters. The second model is based on the PSDFs 
normalized to 20 sec dui’ation and a peak ground acceleration (PGA) value of 1^, 
and does not consider the strong motion duration as parameter. For both models, 
the regression coefficients have been estimated from the PSDFs corresponding to 
the recorded data of western United States for three different geologic site condi- 
tions In case of each recorded time-history, the PSDF is assumed to correspond 
to an ‘equivalent stationary’ ground acceleration process which approximately 
gives the same expected response spectra on considering the transient transfer 
function as those obtained from the time-history The proposed models have 
been illustrated by considering three example ground motions and a paramet- 
ric study has been carried out to discuss the effects of governing parameters 
on PSDF shapes. For crude estimation of normalized PSDFs, this study also 
proposes simple functional relationships for two categories of site conditions. 
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CHAPTER I 
INTRODUCTION 

1,1 General Introduction 

For proper evaluation of the design seismic response of a structure in a 
seismic environment, it is necessary to first characterize the seismic hazard This 
is most commonly accomplished by the specification of (elastic) design spectrum 
curves for different damping ratios. These curves are usually specified m teims 
of an estimated peak ground acceleration (PGA) or effective peak acceleration 
(EPA) at the site corresponding to a ‘scenario earthquake’. Considering the 
seismic environment of the site, this scenario earthqualce is established along a 
nearby fault such that there is a desired level of confidence of the specified ac- 
celeration level not being exceeded in given number of years. The normalized 
design spectrum curves are often obtained after statistical processing of the PGA- 
normalized spectra corresponding to several past ground motions recorded in the 
similar site conditions When large data is available for a seismic environment, 
scaling equations may also be established (e.g , Trifunac and Anderson (1977, 
1978)) and used to directly obtain more realistic estimates of design spectra. 
Since the scenario earthquake approach does not provide a balanced view of the 
seismic hazard when there are several faults near the site tmder consideration, 
balanced risk-based responses of structures have been obtained by integrating 
hazard models for ground motions with the random vibration-based models of 
structural response (see Gupta (1994), Todorovska (1994)). In the scenario earth- 
quake as well as the balanced risk approach, it is desirable to consider the scaling 


2 


of power specUal density function (PSDF) for api)lication to the multi-degrce-of- 
freedom (MDOF) structures. The scaling of acceleration response spectra alone 
as per the conventional practice is undesirable due to the requirement of a modal 
combination rule. Further, by using the response spectra, it is not possible to 
have any information about the higher order peaks, which may play a signifi- 
cant role in the structural damage assessment (Basu and Gupta (1995)), The 
response spectra can also not be used in case of the non-classically damped sys- 
tems due to the fact that the modal response in such a system is strictly not 
analogous to the response of a single-degree-of-freedom (SDOF) system. The 
use of Fourier spectra along with the response spectra accounts only for some of 
these problems A direct scaling of PSDF in terms of various earthquake and site 
parameters, say earthquake magnitude, epicentral distance, focal depth, geologic 
site conditions, etc., appears to be more convenient for the response analyses of 
variety of structural systems. 

To describe a ground motion process through its PSDF, there exist var- 
ious approaches as described by Gupta et al (1994). Analytical descriptions of 
PSDF have been suggested by idealizing the ground motion process to be (a) 
ideal white noise, (b) band limited wliite noise, or (c) filtered white noise pro- 
cess. Ideal white noise idealization is the simplest idealization as it corresponds 
to a PSDF with uniform value for all the frequencies. Band- limited white noise 
idealization offers more realistic improvement over this as here, the PSDF has 
uniform value only within a particular frequency range and is zero outside this 
range. The filtered white noise idealization corresponds to the PSDF proposed by 
Kanai (1957) and Tajimi (1960). This form of PSDF known as the Kanai-Tajimi 
model has been most popular with the researchers for the analytical studies. 
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In this model, the ground acceleration is visualized as the absolute acceleration 
lesponse of a SDOF oscillator subjected to ideal white noise excitation at its 
base. The major advantage in this model is that a wide range of PSDFs can 
be described by simply varying the two parameters, i.e the ground frequency 
and the ground damping. Many researchers (e.g, Lai (1982), Elghadamasi et 
al. (1988)) have suggested the suitable values of these parameters for different 
soil conditions. The main drawback of this model however is that it yields a non- 
zero value at zero frequency which is incompatible with the finite energy m the 
ground motion To correct this, Clough and Penzien (1993) proposed the use of 
one more filter on the Kanai-Tajimi PSDF For the practical applications, how- 
ever, the Fourier spectrum-based PSDFs or the response spectrum-compatible 
PSDFs have been more popular. Even though the earthquake ground motion 
is basically a nonstationary process, its strong motion component may often be 
treated as a weakly stationary process (Elghadamsi et al. (1988)), and the (tem- 
poral) PSDF can then be described in terms of its Fourier spectrum and duration 
(Gupta et al (1994)). For the motions where the initial building up portion and 
the decaying portion at the tail are very short compared to the total duration (as 
in the motions of long duration), this method yields a satisfactory description 
of the ground motion PSDF. As the design ground motions are commonly spec- 
ified through the design (response) spectra, the response spectrum- compatible 
PSDFs (e.g, see Kaul (1978), Sundararajan (1980), Unruh and Kana (1981), 
Pfaffinger (1983) Gupta and Joshi (1993)) are often used in the practical PSDF- 
based applications. The concept of response spectrum-compatible PSDF is based 
on the idea of computing the probabilistic response spectrum from the ground 
motion PSDF (Udwadia and 'IVifanac (1974)), and on modelling the input ac- 



4 


celeraiion process to be an ‘equivalent stationary process’ Such a PSDF is 
computed by an iterative scheme, so as to match the expected peak i espouses of 
a set of SDOF oscillators with the given response spectrum ordinates. In various 
schemes proposed on this basis, it has been found that the PSDF curves obtained 
from response spectra for different damping ratios are substantially different from 
each other. This is due to the basic assumption in these methods that the re- 
sponse of a dynamical system to a stationary process is a stationary process. 
However, it has been shown by Caughey and Stumpf (1961) that, apart from 
the inherent nonstationarity in the ground motion, additional nonstationarity is 
introduced into the response process due to the finite operating time of excita- 
tion. The method proposed by Shrikhande and Gupta (1996) explicitly takes 
care of this additional nonstationarity by the use of transient transfer function. 
In this method, the earthquake ground motion process is modelled as a fiiiite 
duration segment of a stationary random process such that this is equivalent to 
the given nonstationary process in terms of the largest peak responses of a set of 
SDOF oscillators This method is more realistic in the sense that it is indepen- 
dent of the system characteristics and thus similar PSDFs are obtained for the 
spectra for different damping ratios yielding a more consistent description of the 
ground motion. Evolutionary PSDF corresponding to an amplitude-modulated 
stationary process (e.g., see Shinozuka (1970), Liu (1970), Lin and Yong (1987)) 
explicitly accounts for the nonstationarity in the ground motion process through 
the use of a time-dependent envelope function Based on this model, Spanos 
and Vargas Loli (1985) have suggested the computation of spectrum- compatible 
evolutionary PSDF. However, this model is not accurate enough due to ignoring 
the frequency nonstationarity in the ground motion, and is not as simple and 
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elegant as the model of equivalent stationary process 

Except for the availability of a few empirical forms of design PSDF in 
the recent years (Chang et al (1991)), no attempt has been made yet to es- 
tablish an empirical relation to scale the ground PSDF m terms of the earth- 
quake parameters hke epicentral distance, magmtude, geologic site condition, 
focal depth etc A substantial amount of work has however been done to de- 
velop scaling relationships for other functionals like PGA, Fourier spectrum, 
pseudo spectral velocity (PSV), and spectral acceleration (SA). Scaling rela- 
tions for PGA in terms of magnitude and epicentral distance have been pro- 
posed by several researchers (for example, see Gutenberg and Richter (1942, 
1956), Neumann (1954), Housner (1965), Blume (1965), Kanai (1966), Milne 
and Davenport (1969), Esteva (1970), Cloud and Perez (1971), Schnabel and 
Seed (1973)), Each of these correlations was obtained by fitting an attenuation 
curve to a limited number of data points. Later, Trifnnac and Brady (1976) 
used the Richter attenuation function (Richter (1958)) to correlate the means 
and standard deviations of the peaks of ground acceleration, velocity and dis- 
placement With magmtude for different site conditions The Richter attenuation 
function was obtained empirically from the data for Southern California, and 
this thus contained an ‘average’ information about the properties of wave prop- 
agation regarding geometric spreading, scattering, etc through the crust in 
Southern California Trifunac (1976a) later improved this model by attaching 
appropriate confidence levels to the predicted peak amplitudes and by consider- 
ing magmtude, site condition, epicentral distance as the scaling parameters in 
the equation By obtaining independent estimates (based on source mechanism 
theory) of the peaks of acceleration, velocity and displacement and by comparing 
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those With the estimates of the proposed scaling relationship, Trifnnac (1976a) 
observed that the attenuation law of Trifunac and Brady (1976) could be ex- 
trapolated to give satisfactory results for distances less than 20 km also The 
model of Trifunac (1976a) was later generalized to other spectral quantities, e g , 
to Fourier spectra by Trifunac (1976b), to absolute acceleration spectra by Tri- 
funac and Anderson (1977), and to pseudo relative velocity spectra by Trifunac 
and Anderson (1978) at a set of discrete periods These models facilitated the 
direct description of design ground motion in terms of the known parameters at 
a desired confidence level. 

The attenuation function used by Trifunac and Brady (1976) had the lim- 
itation that it was independent of parameters hke magnitude (which effectively 
represented the source dimension of an earthquake), focal depth and geological 
environment of the site. In the early 80’s, as more strong motion data from 
Northern and Southern Cahforma became available, it was attempted to model 
the attenuation in a better way by incorporating the fact that the high-frequency 
waves attenuate much faster compared to the low frequency waves. Trifunac and 
Lee (1985a) suggested a frequency-dependent attenuation function which was 
composed of two parts, one depending on the frequency, i.e the attenuation 
coefficient, and the other depending on the source- to-station representative dis- 
tance. The ‘representative distance’ as proposed by Gusev (1983) varies as a 
function of the earthquake magnitude, and thus incorporating this in the atten- 
uation function represented an improvement over the Richter model. The new 
attenuation function was first employed for the scahng of Fourier amplitude spec- 
trum by Trifunac and Lee (1985b), and later in several other scaling studies (e.g , 
see Trifunac and Lee (1985c), Trifunac (1987), Lee (1987)). 
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In various scaling equations presented in the past 20 years (e g , see Tn- 
funac (1976a, 1976b), Trifunac and Anderson (1977), Trifunac and Lee (1978, 
1979, 1985b, 1985c)), the rate of growth of spectral amplitudes with earthquake 
magmtude has been considered to be linear up to a certain minimum value of the 
magmtude Beyond this range and up to a certain maximum value, the growth is 
assumed to have a parabohc form The spectral amplitudes are considered to be 
invariant of the magmtudes higher than this maximum value. Though this func- 
tional form of the growth lacked any physical justification, it was supported by 
the observed trend of data (Trifunac (1976a, 1976b), Trifunac and Brady (1976)). 

It IS well known that the various ground motion functionals have strong 
dependence on the geologic and soil conditions at the site. In some of the works 
reported (eg., those by Trifunac (1976b), IVifunac and Anderson (1977), Tri- 
funac and Lee (1985b, 1985c)), the site classification was based on the geological 
characteristics in terms of the hardness of the soil surrounding the recording sta- 
tions This classification was proposed by Trifunac and Brady (1976), and the 
main advantage of using such a classification is that this information is available 
from the knowledge of surface geology only. However, in this classification, the 
‘geometric characteristics’ (i.e., the horizontal and vertical extent of the local 
inborn ogenei ties) regarding the local soil conditions are overlooked To overcome 
this drawback, the recording site condition was described through the depth of 
sediments beneath the recording station (Trifunac and Lee (1978)). The idea be- 
hind this basis of site classification was that the wave velocity, material rigidity 
and density increase substantially in going with depth from sediments to rock 
Later, in a further improvement, the parameter of local soil conditions as de- 
fined by Seed et al (1976) was considered simultaneously with the geological site 
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conditions by Tnfnnac (1987) and Lee (1987) 

It IS apparent from the above that more and more parameters were in- 
cluded in the scaling relations with the passage of time for narrowing down the 
scatter of data points. Further, in this direction, Lee et al. (1995) has recently 
presented scaling relations for peak ground acceleration, velocity and displace- 
ment by explicitly including the effect of the transmission path 

Along with the magmtude-based scaling relationships as discussed above, 
relationships incorporating earthquake intensity, e g , Modified Mercalli Intensity 
(MMI), as a parameter have also been developed for different spectral quan- 
tities In spite of having inherent drawbacks, e.g , large uncertainty in the 
intensity data, these relationships have been developed because in many seis- 
mic regions of the world, no measurements of the strong ground motions are 
available Many historical records are in fact based on locally developed inten- 
sity scales only The intensity-based scaling relationships were first proposed 
for scaling peak ground acceleration (for example, see Ishimoto (1932), Kawa- 
sumi (1951), Hershberger (1956), Richter (1958), Neumann (1954), Medvedev 
and Sponheuer (1969), Okamoto (1973), Savarensky and Kirnos (1955), Trifunac 
and Brady (1975a)). Later, scaling relationships were also developed for other 
functionals, e.g , by Trifunac (1979) for absolute acceleration spectra, by Tri- 
funac and Anderson (1978), Trifunac and Lee (1985b), and Trifunac (1987) for 
Fourier spectrum, and by 'IVifunac and Lee (1978), Trifunac and Lee (1985c), 
and Lee (1987) for PSV. 

As mentioned earlier, more refined scaling relationships with greater 
number of parameters have been able to describe the ground motion functionals 
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more accurately as has been indicated by the gradual narrowing down of the 
spread of the spectral amplitude over a given confidence interval It is expected 
that with further increase in the size of the database, further refinements will 
be proposed in the scaling models involving perhaps more number of parameters 
like transmission path, stress drop, seismic moment, etc. However, due to the 
practical constraints, all of these parameters may not always be available to the 
designer in seismic hazard quantification Keeping this in view and for simplicity 
in the first step towards developing a scaling model of PSDF, the scaling model 
as proposed by Trifunac and Lee (1985b) for Fourier spectrum amplitudes has 
been considered as the basis in this study This model makes use of those pa- 
rameters which are routinely employed to describe the overall size and seventy 
of the earthquake, namely, epicentral distance, magmtude, focal depth, and ge- 
ologic site condition Further a model for the Fourier spectrum amplitudes has 
been considered as the basis due to the direct dependence of temporal PSDF of 
a ground motion record on its Fourier spectrum amplitudes and strong motion 
duration (see, for example, Bendat and Piersol (1986)) 

Two scahng models have been presented in this study. Both models con- 
sider epicentral distance, magnitude, focal depth, and site geology as the govern- 
ing parameters. In addition to these, the strong motion duration has also been 
considered as a parameter in case of the first model. Since there is no standard 
and well accepted definition of strong ground duration, though several definitions 
have been proposed in the past two decades, e g , those by Bolt (1973), Page et 
al (1972), Trifunac and Brady (1975b), McCann and Shah (1979), Kawashima 
and Aizawa (1989), Vanmarcke and Lai (1980), Mohraz and Peng (1989), and 
Trifunac and Westermo (1982), the second model has been proposed for ‘normal- 
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ized’ PSDF to a fixed duration of 20 sec The PSDF data used for the second 
model IS also normalized to a fixed value of expected PGA, besides normalization 
to the same duration, since it is common to characterize the severity of design 
ground motions by PGA values The PSDF data has been generated from the re- 
sponse spectra of the recorded motions by using the spectrum-compatible PSDF 
definition of Shrikhande and Gupta (1996) Both of the proposed models do not 
consider the site geology as a parameter Instead, the scaling coefficients have 
been regressed separately for three different (geologic) site conditions Also, both 
models have been proposed for the horizontal component of ground motion only, 
and thus, the parameter of component direction also has not been considered. 
With a view to further simplify the model for normalized PSDF for crude ap- 
plications, simple average shapes have also been estimated for normalized PSDF 
on the same hnes as in Seed et al. (1976) for PGA-normalized response spectra 
The database used for the present study consists of the 494 ‘free-field’ records 
of 106 earthquakes in Western United States, starting from 1931 Long Beach 
Earthquake, California to the 1984 Morgan Hill Earthquake, Cahfornia, as in 
Lee and Trifunac (1987) 

1*2 Organization 

This work has been presented in three chapters following this chapter 

In Chapter II, the first scaling model, correlating the horizontal ground 
PSDF with earthquake magmtude, epicentral distance, focal depth and duration, 
has been presented for different site conditions. Least square estimates for vari- 
ous period-dependent coefficients have been obtained by regression analysis and 
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the deviations of the actual PSDF values from the estimated values have been 
described probabilistically The proposed model has been illustrated by com- 
paring the PSDFs of a few example recorded motions with their probabilistic 
estimates as obtained from the proposed model, and by presenting the estimated 
PSDFs with 50% level of confidence for varying durations, magnitudes, epicentral 
distances and site conditions 

In Chapter III, the second scaling relationship and the ‘average’ PSDFs 
based on the ‘normalized’ PSDF data have been presented As m Chapter II, the 
second scaling relationship has also been illustrated through comparison of actual 
PSDFs with the probabilistic estimates of PSDFs and through a parametric 
study Further, simple functional forms have been suggested for the average 
PSDFs m different site conditions 

A brief summary and conclusions of this study have been presented in 


Chapter IV. 
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CHAPTER \1 

SCALING OF PSDF IN TERMS OF EARTHQUAKE 
MAGNITUDE, EPICENTRAL DISTANCE, FOCAL 
DEPTH AND STRONG MOTION DURATION 


2.1 Brief Overview 

In this chapter, a new regression model has been proposed for scaling the 
PSDF of the horizontal component of ground motion process. This assumes the 
parametric dependence of PSDF on earthquake magnitude, epi central distance, 
focal depth, strong motion duration and site geology. The proposed model is 
a simple extension of the model proposed by Trifunac and Lee (1985b) for the 
scaling of Fourier spectrum amplitudes. The PSDF corresponding to a recorded 
ground motion has been estimated from the response spectra of this motion by 
using the procedure suggested by Shrikhande and Gupta (1996) The scaling 
coefficients have been determined from the PSDFs corresponding to a data-set 
of 988 recorded ground motions for three different geological site conditions. 
These ground motions have been recorded during the 106 earthquake events in 
the Western U.S.A region till 1984. In order to use the regressed coefficients 
for estimating PSDF at a given confidence level, the inherent variations in the 
recorded PSDF values have been probabilistically modelled at different periods. 

2.2 Scaling Relationship 


In. the absence of any previous study on the scaling of PSDFs, a scaling 
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ieIa(ioiis[iip consistent with the following relationslup between Fourier spectrum, 
FiS'(T), and PSDF, PSDF(T)y has been considered According to tliis, the 
(temporal) PSDF of the ground motion at period, T, may be estimated as (e g., 
see Bendat and Piersol (1986)) 

PSDF{T) = . (2.1) 

TT jTfl 

where, Tg is the strong motion duration of the ground motion. On considering 
the preliminary scaling model for Fourier spectrum as presented by 'Ildfunac and 
Lee (1985b), the scaling equation for PSDF of the horizontal ground motion may 
be assumed to be 


logio PSDF(T) = 2M + 2An{A,M,T) + bi{T)M + 2>2(T)A + b 3 (T) 

+ b,(T)M^ + h{T) log, aT,. (2.2) 


More sophisticated models for scaling of the Fourier spectrum amplitudes also 
consider a parameter for the local site conditions (Trifunac (1987)). However, 
this model has been chosen in order to obtain a preliminary scaling relationship 
which uses only the readily available earthquake paramaters. In Eq (2.2), fcj(T')’s 
are the period-dependent coefficients and M is the ‘published’ magnitude. A 
‘published’ magnitude corresponds to the local magnitude if iW <6.5 and to the 
surface wave magnitude if M >6.5 (Trifunac and Lee (1985a)). Further, A is the 
‘representative distance’ which basically represents the effective source-to -station 
distance It has been defined by Gusev (1983) in terms of the epicentral distance, 
i?, the focal depth, the size of the fault, i?, and the coherence radius, of 
the source function as 





Ri + H'^ + Soy 



(2 3) 
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The coherence radius, 5o, is period- dependent and it can be appioxiniated as 
(3Tl2 where P is the velocity of radiation of waves having period, T (Trifunac 
and Lee (1985a)). has been assumed to be 1000 m/sec in this study The term, 
T) is the period-dependent attenuation function proposed by Tidfunac 
and Lee {1985a) as 

=A(T)logioA R<Rq 

= A)(T)logioAo-(i?-i?o)/200 R>Rq, (2.4) 

where is the transition distance from the earthquake source at wliich the 
surface waves start dominating the ground motion, and Aq is the representative 
distance corresponding to Rq FWther, Aq{T) is the empirically determined 
attenuation function at time period, T, and is represented as (Trifunac and Lee 
(1985a)) 

Aq(T) = -0.767079 + 0.271556 logior - 0 52564 (logioT)^ ; T < 1.8 sec 
= -0.732025 ; T > 1.8 sec . 

(2.5) 

As the values of the fault size are usually not available, an empirical relation- 
ship has been proposed by Trifunac and Lee (1985a) to estimate this from the 
earthquake magnitude at each period as 

S' = 0.2 + - 1 (S6.5 -02), (2,6) 

where 56 5 is the fault size for M =6.5. Its value for six different period ranges 
has been obtained by Trifunac and Lee (1985a) through iterative calculations, 
and those values have been adopted for the present study. 
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It. raay be noted that the prc.se nt study considers the horizontal gionnd 
motions only. Hence, the parameter, v which has been used by Trifunac and 
Lee (1985b) for the component direction of motion has been omitted The pro- 
posed model also considers the effect of geologic site conditions in a different 
manner. Instead of considering s as a parameter in regression analysis, with 
s = 0, 1, 2 representing respectively the alluvium, intermediate and hard rock 
conditions, three separate regression analyses have been carried out for these site 
conditions. Thus, instead of single scaling equation with s included as a param- 
eter, three different scaling equations have been developed in this study for the 
tliree different site conditions It raay also be noted that the new scaling param- 
eter, Tgy does not have a unique or well-accepted definition since the earthquake 
ground motions are highly non-stationary and Tg is deemed to refer to the strong 
motion part of an earthquake ground motion. For the present study, the defini- 
tion by Trifunac and Brady (1975b) has been adopted since the strong motion 
duration according to this definition can be derived from a recorded accelerogram 
by considering the arrival of seismic energy with time in a very simple manner. 
Due to this reason, this definition has been widely used in several studies over 
the past two decades. By this definition, the continuous time interval in which 
the most significant contribution to the seismic energy takes place is considered 
to be the ‘strong motion duration’. The seismic energy at a time instant, t = f', 
is measured by the integral of square of acceleration from t = 0 to t = and 
the strong motion duration is defined as the difference of those time instants at 
which this integral attains 5% and 95% of the tgtal seismic energy (measured at 
the end of the record). 
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2.3 Raw PSDF Data and Data Selection 

The present study is based on the EQINFOS data bank (Lee and Tri- 
funac (1987)), This consists of 988 records for horizontal component of ground 
motion from 106 earthquakes which were recorded from 1931 Long Beach Earth- 
quake, California to the 1984 Morgan Hill Earthquake, in the Western U.S.A 
region Some of these records have been left out of the present regression anal- 
ysis due to the lack of complete details. For example, no reliable estimates of 
earthquake magnitude were available for AU294, AU295, AU297, AU298 records. 
For the event ^ 7, 10, 11, 12, 25, 34, 35, 38 and for record ^ 1341 and 1381, the 
recorded time histories were incomplete. For some of the records, namely BA142, 
AB025, AU295, AU297, AB032, AB029, AB028, AU310, BA144, BA255, AY372, 
no information about the geologic environment of the site was available. Thus, 
this study has been based on 960 recorded accelerograms, 733 for alluvium site 
conditions, 85 for intermediate site conditions and 142 for hard rock site con- 
ditions For a few early records, the focal depth, R, values of the earthquakes 
were not available. In such ceises, H has been uniformly assumed to be 16.0 km 
(Trifunac and Brady (1975a)). 

For evaluating PSDF corresponding to each accelerogram record prior 
to the regression analysis, the method of computing spectrum- compatible PSDF 
as proposed by Shrikhande and Gupta (1996) has been employed with the du- 
ration taken to be the same as the strong motion duration of the record. Three 
smoothed response spectrum curves corresponding to 0,01, 0.02 and 0.05 damp- 
ing ratios have been considered in case of each record, and these spectra have 
been assumed to be the ‘expected’ spectra for the underlying ground acceleration 
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pi o cess 111 each of Ihc three cpscs, the matching of the computed spectrum with 
the target spectrum has been afhieved within 1 5% of the average error at 91 
periods between 0.04 and 15,0 sec. The three PSDFs obtained from the tliree re- 
sponse spectra have been then averaged out to give the mean information about 
the PSDF of the underlying process. Let tliis be denoted by PSDF{T). 

Since the number of records for different magnitude classes, namely 3 0- 
3.9, 4 0-4,9,. 7. 0-7.9, were substantially different from one another, and for 

some of the earthquakes (for example, 1971 San Fernando Earthquake), there 
were excessive number of records which could have a bias on the final coefficient 
values, data screening has been earned out. To accomplish this, a technique as 
reported in the previous works on regression analyses (e g , see Trifunac (1976b), 
Tnfunac and Anderson (1977, 1978), Trifunac and Lee (1978, 1985b, 1985c)) 
has been employed. For each of the site categories (i.e , alluvium, intermedi- 
ate and hard rock site), all the PSDF data has been partitioned into 5 groups 
corresponding to the different magnitude classes. Further, to properly balance 
the effects of attenuation at small and large distances, the data in each of the 
above groups has been sub-divided into two sub-groups, one for 72 < 100 km 
and the other for 72 > 100 km. Various PSDF data points in each of these 
sub-groups have been then rearranged at a period, T, in such a fashion that the 
value, {\ogiQPSDP(T) — 2M — 2Att{A^M^T)) decreased monotonically. Since 
the objective of the data selection procedure is to ensure that approximately 
same number of data points represent each of the sub-groups, different values 
of confidence level, say p, have been chosen and the data points in a sub-group 
which correspond to those levels only have been considered Thus, if m levels of 
probability, p = Pi, P 2 )'-) Pmj are specified and if fci, fc2v”i denote the nearest 
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integers to pin^ p^n respectively, the A:ith, A:2th, k^th data points in a 

sub-group of n data points have been considered Thus, in case of the sub-groups 
coriesponding to < 100 km, the chosen p values are pi = 0 05, P2 = 0.10, , 
pm = 0 95 With m = 19 In some of these sub-groups, if n is less than 19, all the 
data points have been taken into account For the sub-groups corresponding to 
R > 100 km, fewer records were available, and hence, a maximum of 5 points 
only have been selected in each sub-group (with pi = 0 1667, p2 = 0 3333,..., ps 
= 0 8333) It may be mentioned here that the above selection process screens 
out different sets of data points at different periods, and thus, all the available 
PSDF data is effectively utihzed in the regression 

2,4 Regression Analysis and Distribution of PSDF 

Linear regression analyses have been carried out on 733 PSDFs for al- 
luvium site conditions, 85 PSDFs for intermediate site conditions and on 142 
PSDFs for hard rock sites at 91 values of T, ranging from 0 04 to 15.0 sec. In 
case of each of these site conditions, let the least square estimates of the coef- 
ficients, bi(T), b 2 (T), b 3 (T)j b 4 (T) and 65 (ff*), as obtained froni the regression 
analysis at each period, T, be denoted by bi(T), b 2 (T), b 3 (T)^ b 4 (T) and &6(T) 
respectively. The 95% confidence intervals of these coefficients have been esti- 
mated by assuming those to be distributed as per the student’s t-distribution 
with (n — 5) degrees of freedom where n is the total number of data points used 
for regression. The variance of each of the coefficients has been computed for 
this purpose as in Westermo and 'Ibifunac (1978). The least square as well as 
the 95% confidence estimates for each coefficient have been smoothed along the 
logioT axis 
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t' 


On feubsLituting the (smoothed) values of the coefficients, t’s, in Eq. (2.2), 
we obtain the least square estimate of PSDF{T) as 


\og^Q FSDF(T) = 2M -b2v4tt(A,M,T) + &i(T)M + b2(T)A + kiT) 

+ h{T)M^ + k(T)\ogiaT,. (2.7) 


For a given set of values of T, Af , A and Tg, logig PSDF{T) represents a parabolic 
variation in M This parabolic dependence on M has been however found to be 
valid, m case of Fourier spectrum amplitudes, only within a particular interval 
of magnitude (Trifunac (1976a)). In fact, for the small magnitudes, the growth 
of logarithm of Fourier spectrum has been found to be proportional to the mag- 
nitude, and for higher values of M, this has been found to be invariant of M 
(Trifunac (1976a)) To model this behaviour in a simple manner, it has been 
assumed in view of Eq (2 1) that the growth rate of logic -E‘S'iPE(T’) ^ 
corresponds to a slope equal to 2 for Af < For Mni\n < ^ 

slope decreases gradually and finally becomes zero for M > be 

noted that this assumption is based on neglecting the variations in strong motion 
duration with magnitude in these ranges. Thus, Eq. (2.7) has been considered 
to be valid only in the range, Mmjn ^ Af < where 


•^mln — 


k{T) 
A(t) ' 


£LIld 


2+k{T) 

2h(T) 


( 2 . 8 ) 


(2.9) 
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Eq. (2 7) is then modified to 

log,oPS5i5’(r) = 2^tt(A, Af,T)+ 

' 2M + k{T)M^,, + S2(r)A + kiT) 

+b4{T)Ml,„ + k{T) logio T, ; M< M„,„ 

^ ■2M + k{T)M + k{T)A + h{T) ^ loi 

+k{T)M^ + k{T)[og,oT,, 

2M,„^ + k{T)M^^ + h2{T)A + k(r) 

, +b4(T)M^^ -h bs(T)logji)Tg ; Mnax < A7 

and the corresponding residual describing the deviation of the recorded PSDF 
from the estimated PSDF at period, T, is given by 

e(T):=logioPSDF(T)~logioJ^sBF(T) . (2.11) 


At a fixed T, all the residuals have been computed from the raw PSDF data 
for each of the soil conditions by using this equation, and the actual probability 
distribution, p*(e, T) has been obtained from the fraction of residuals. According 
to this, the piobability that a chosen value of 6(T) will not be exceeded has been 
evaluated simply by finding out the ratio of the number of residual values less 
than the chosen e(T) to the total number of residual data values. It is assumed 
that a normal distribution function with the mean, p(T), standard deviation, 
o-(T), and described as 


p(e,r) = 



exp 



( 2 . 12 ) 


closely approximates the computed distribution, p*(e,T). Hence, the values of 
mean and standard deviation, say /i(T) and d(T) respectively, which provide 
the best fit of Eq, (2.12) withp*(e,T') have been estimated, These least square 
estimates, /i(T) and S‘(T), have been then smoothed along logjoT axis for each of 
the soil conditions. The smoothed values, /i(T) and d(r), now form the basis for 
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estimaliiig the PSDF with a desired level of confidence. At a period T, Eq (2.12) 
can be used to estimate e{T) corresponding to a given level of confidence and be 
added to the calculated PSDF{T) (see Eq (2 10)) for this purpose. To check 
the reliability of the assumption that the residuals are normally distributed, two 
statistical ‘goodness of fit’ tests, namely the chi-square test and the Kolmogorov- 
Smirnov (KS) test, have been conducted. The KS statistics, KS{T) and the chi- 
square statistics, X^(^) have been calculated at each of the 91 periods for different 
site conditions, and have been compared with the respective 95% cut-oflP values. 

2.5 Results and Discussion 

The smoothed least square estimates, hi{T) through bs{T)y and the 95% 
confidence estimates of the regression coefficients are shown in Figs. 2, 1(a) to 
2 1(c) for alluvium site (s = 0), intermediate site (s = 1) and hard rock site (s 
= 2 ) conditions respectively. The least square estimates are shown by the sohd 
lines while the 95% confidence estimates are shown by the dashed lines. Tables 
2.1(a) to 2 1(c) present the (smoothed) least square estimates, 6i(T) through 
bs{T)y and the (smoothed) values of /i(T), c{T) at 12 selected periods for the 
three site conditions. 

Figs. 2.2(a) to 2 2(c) show the comparison of the actual residuals with 
the estimated residuals (corresponding to the smoothed fi{T) and d-(!r) values 
in Eq. (2,12)) for the alluvium site, intermediate site, and hard rock conditions 
respectively. In each figui'e, the solid lines represent the actual residuals for 
p* = 0.1, 0.3, 0.5, 0.7 and 0.9 while the dashed lines represent the estimated 
residuals for p = 0.1, 0.3, 0.6, 0,7 and 0.9. It may be observed that the estimated 
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Table 2 1(a) - I^nrisl Stjuriu' Esliinatc's of Rcgiession (’offlit u-iils and Rr'iidua] 
Pfiianiolcis fo] Alluvium SilP Condition 


Pciiod, 

Least Square Estimates of Regression Coefficients 

Ro.sidual Paiametci's 

T (see) 

h(T) 

10062(T) 

h(T) 

JO^(T) 

h(T) 

KT) 

om 

0 04 

-1 720 

1.480 

I 910 

0.735 

-0.8II 

0 109 

0.673 

0 065 

-1 830 

1.100 

2 220 

0 757 

-1 050 

0.104 

0 650 

0 11 

-1 840 

0 725 

2 150 

0 749 

-1.130 

0.109 

0 640 

0.19 

-1.750 

0 435 

1 370 

0 724 

-0.974 

0 117 

0 638 

0.34 

-1 690 

0 274 

0 224 

0 762 

-0.653 

0.120 

0 663 

0.50 

-1 680 

0 175 

-0 590 

j 0 797 

-0 353 

0 121 

0.707 

0.90 

-1 680 

-0 175 

-2.020 

0 871 

0 376 

0 136 

0 822 

1 60 

-1.810 

-0.560 

-3.330 

1 060 

1 260 

0.167 

0 979 

2.80 

-2 080 

-0.464 

-4.260 

1.340 

1.940 

0.199 

1,213 

4.40 

-2 330 

0 013 

-4 710 

1.540 

2.250 

0.189 

1 443 

7 50 

-2.540 

0,444 

-4.700 

1 540 

2.230 

0 075 

1,583 

14 00 

-2 640 

0 379 

-4 530 

1 360 

2 040 

-0.101 

1.517 











26 


71ahle 2.1(b) - boasl Squaio of Hogiossjon Coofhcn'iith anci r?osidiia] 

Paiameters foi Intoi mediate Site Condition 


Pciiod, 

1 

Lca.st S( 

]uaie Estiiiia 

tes of Pegie.ssion Coefliucnls 

Residual PaianicLeis 

T (sec) 

bi(T) 

1 lOO^zCT) 

h{T) 

m^iT) 

k{T) 

/i(T) 

HT) 

0.04 

0 179 

0 590 

-3 050 

-0 832 

-0 971 

-0 015 

0 395 

0 065 

0 257 

0 442 

-3 130 

1 -0 991 

- 1.430 

-0 007 

0 418 

1 

Oil 

0.575 

0.156 

-3 960 

-1 310 

-1 650 

-0 007 

0 437 

0 19 

0 845 

-0 280 

-5 130 

- 1.510 

-1 500 

- 0.011 

0.441 

0.34 

0 708 

-0 759 

-5 710 

- 1.280 

- 1.130 

- 0.011 

0 450 

0 50 

0 492 

-1 060 

-5 900 

-1 000 

- 0.856 

- 0.013 

0 477 

0.90 

0 078 

- 1.510 

-6 310 

- 0.482 

- 0.312 

- 0.020 

0 564 

1,60 

-0 923 

- 1.800 

- 5.550 

0.534 

0 541 

- 0.021 

0.699 

2 80 

-2 310 

- 2.480 

-3 650 

1.860 

1 430 

- 0.020 

0.902 

4.40 

-3 420 

- 2.800 

-1 780 

2,820 1 

1 950 

- 0.040 

1.094 

1 7.50 

- 4.690 

- 1.940 

1 no 

3 640 

2.060 

-0 129 

1.257 

14.00 

- 5.710 

-0 529 

3.930 

4.080 

1 790 

- 0.302 

1.291 
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Table 2 1(c) - Least Sciuaie Hstiindlch of Rogiossioii Coolfkionls and Residual 
Para ni del’s for Raid Rock Site Condition 

Pciiod, Least Squaie Estimates of Regicssion CoefRuents Residual Paiametcrs 
r(sec) hi{T) mker) b,{T) \ok{T) k{T) A(r) a{T) 

0 04 -1 850 1 730 2 030 0 911 -0.836 -0 015 0 651 

0 065 -1.200 1.410 0 526 0 233 -1 240 -0.042 0 664 

011 -0 882 1 000 -0 540 -0 073 -1 210 -0 051 0.688 

0 19 -0.398 0 622 -2 530 -0 527 -0 682 -0.055 0.750 

0.34 0 550 0.000 -6 110 -1 450 0 298 -0.063 0.859 

0 50 0,718 -0.696 -7 450 -1 560 0.898 -0 071 0.957 

0.90 0 006 -1.600 -7 210 -0.785 1 840 -0.089 1.148 

1 60 -0 991 -1 510 -6 230 0.192 2 600 -0.109 1 342 

2.80 -1,900 -1.100 -4 930 1 140 2 720 -0146 1 509 

4,40 - 2.630 -0.961 -3.840 1 810 2 420 -0 208 1 615 

7.50 -2 780 -1.150 -3.560 1.840 1 790 -0.356 1.666 

14.00 -2 710 -1 220 -3.550 1.590 0 960 -0.572 1.706 











Residues e(T) Residues e(T) 



Figure 2.2(a) Residuals forp,p* = 0.1, 0.3, 0.5, 0.7 and 
0.9 in Case of Alluvium Site Condition. 



Figure 2.2(b) Residuals for p^ p* = 0.1, 0.3, 0.5, 0.7 and 

0.9 in Case of Intermediate Site Condi- 
tion. 





Actual 

Estimated 


Time Period, T (sec) 

Figure 2.2(c) Residuals values for p, p* = 0.1, 0.3, 0.5, 
0.7 and 0.9 in Case of Hard Rock Site Con- 




30 


residuals represent nice and smooth approximations for the actual icsiduals. The 
estimated residual values for nine different probability levels have been presented 
in Tables 2 2(a) to 2 2(c) respectively for the three site conditions 

The KS statistics, KS{T)y and the chi-square statistics, X^(^)) ^^-ch 

period, T, have been compared with their respective 95% cut-off levels (repre- 
sented by the dashed lines) for all the three site conditions in Figs. 2 3(a) and 
2 3(b) respectively. From Fig. 2.3(a), it is seen that except for the long periods, 
the computed KS statistics are well within the cut-off level for the alluvium site 
condition case. The KS test has however passed at fewer points for the other two 
site conditions. Further, Fig. 2 3(b) shows that for intermediate site conditions, 
the results of chi-square test are in very good agreement with the hypothesis. 
For the hard rock case also, a good agreement exists except in case of a few high 
periods. The chi-square test has however failed to accept the hypothesis to be 
true in case of alluvium site conditions. It is seen from these results that the nor- 

i 

mat distribution may not be an ideal description of the error distribution at all 
time periods. However, since the hypothesis has peissed through both the tests 
at several time periods, and since there appears to be a very good agreement 
between the actual and estimated residuals (see Figs. 2.2(a) to 2 2(c)), the error 
estimates based on the assumed normal distribution have been considered to be 
acceptable. 

In order to compare the probabilistic estimates of PSDF with the actual 
PSDFs corresponding to the recorded data, three example cases have been con- 
sidered. These correspond to (a) east component recorded at Gilroy 3 station 
during the 1984 Morgan Hill Earthquake with M = 6 2, i? = 38,9 km, H — 9.0 



Table 2.2(a) - Residual Values for Alluvium Site Condition 
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Table 2.2(b) - Residual Values for Intermediate Site Condition 
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Table 2.2(c) - Residual Values for Hard Rock Site Condition 
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Figure 2.3(a) Comparison of KS Statistics at Different 
Time Periods with the 95% Cut-off Level 
for Different Site Categories. 
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Figure 2.3(b) Comparison of Chi-square Statistics at 
Different Time Periods with the 95% Cut- 
off Level for Different Site Categories. 
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km, 5 = 0 and = 21 08 sec, (b) south component recorded at Southern Cali- 
fornia Edison Company, Colton station during the 1970 Lytle Ci'cek Earthquake 
with M ~ 5 4, i? = 31 5 km, H — SO km, s = 1 and Tg = 10.60 sec, and (c) 
SC5E component recorded at Wrightwood station during the 1971 San Fernando 
Earthquake with M — 6.4, i? = 70 6 km, H = 13 0 km, s = 2 and Tg = 15 98 sec. 
Figs, 2 4(a) to 2.4(c) respectively show the comparisons for these cases, with solid 
lines representing the estimated PSDFs corresponding to the probability levels, p 
= 0.1, 0,5 and 0.9 in each figure, and dashed lines representing the actual PSDFs. 
The estimated PSDFs have been obtained by using the least square estimates 
and error estimates as given in Tables 2 1(a) to 2 1(c) and Tables 2.2(a) to 2.2(c) 
respectively at 12 time periods. The region bounded by two extreme solid lines 
in each of the three figures represents the respective 80% confidence interval. It 
is seen that the actual PSDFs he well within these confidence intervals Also, 
in the second and third example cases, the actual PSDFs are in good agreement 
with the estimated PSDFs for p = 0.5. However, the large uncertainty in the 
estimation of PSDFs as implied by the observed wide confidence intervals indi- 
cates that the proposed model needs to be improved further possibly with the 
inclusion of more earthquake source, path and site condition parameters. 

To illustrate how PSDF changes with the variations in strong motion 
duration, earthquake magnitude, epicentral distance, and geologic site condition, 
Figs. 2,5 to 2.8 have been obtained by using the least square estimates and 
error estimates as given in Tables 2, 1(a) to 2,1 (c) and Tables 2.2(a) to 2,2 (c) 
respectively at 12 time periods, and by taking focal depth, H — ^ km and p 
= 0.5, Fig. 2.5 shows the PSDFs for Tg = 10, 15, 20 and 25 sec in case of 
alluvium site conditions, with M and R respectively kept constant at 6 5 and 




Figure 2.4(a) Comparison of Actual and Estimated 

PSDFs for Morgan Hill Earthquake 
Case. 



. Figure 2.4(b) Comparison of Actual and Estimated 

PSDFs for Lytle Creek Earthquake 
Case. 





Figure 2.4(c) Comparison of Actual and Estimated 
PSDFs for San Fernando Earthquake 
Case. 
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Figure 2.5 Estimated PSDFs for Different Durations 
in Case of Alluvium Site with R = 50.0 
km, if = 5.0 km, Tg = 20.0 sec and p = 





Estimated PSDF ( 


O 2 

0 ) 10 ® 


10 


0.1 

Time Period, T (sec) 

Figure 2.6 Estimated PSDFs for Different Magni- 
tudes in Case of Alluvium Site with. R = 
50.0 km, = 5.0 km, Tg = 20.0 sec and 
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Figure 2.7 Estimated PSDFs for Different Epicentral 
Distances in Case of Alluvium Site with 
M = 6,5, H = 5.0 km, Tg = 20.0 sec and 
p = 0.5. 
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Figure 2*8 Estimated PSDFs for Different Site Con 
ditions with M = 6,5, = 50.0 km, H = 

5.0 km, Ts == 20.0 sec and p == 0.5. 
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50 km Fig 2 C shows the PSDFs for M = 4 5, 5.5, 6.5 and 7 5 in case of 
alluvium .site conditions, with R and Tg respectively kept constant at 50 km 
and 20 sec. Fig 2 7 shows the PSDFs for R = 25, 50, 100 and 200 km in 
case of alluvium site conditions, with M and Tg kept constant at 6 5 and 20 
sec Similarly, Fig 2 8 shows the PSDFs for the three site conditions, with the 
other parameters taken asjy^ = 65,ii = 50 km, and Tg ~ 20 sec. In any of 
these figures, iiowever, a mutual comparison of various curves docs not give a 
strictly true picture due to the observed correlations of strong motion duration 
with the other earthquake and site parameters (see iVifunac and Brady (1975b), 
Westermo and Trifunac (1978), Novikova and Trifunac (1994)). For example, 
in case of Fig. 2 6, it is unlikely that the events with different magnitudes will 
lead to 20 sec of strong motion duration with no change in the other parameters. 
In fact, it has been observed that the large magnitude earthquakes may cause 
longer ground motions at a site and thus, a comparison of various curves in 
Fig. 2 6 will be complete only when a higher value of Ta is associated with a 
larger value of M. Observing from Fig. 2.5 that the PSDF for a longer strong 
motion duration ground motion is likely to be richer in longer period waves, the 
shifting of the peaks to the longer periods with increasing M as in Fig. 2.6 may 
be actually more pronounced. This is consistent with the observation of Trifunac 
and Lee (1985b) that large magnitude earthquakes are dominated by the longer 
period waves. Same trend should also be obtained with the increasing epicentral 
distances from a simultaneous look at Figs. 2.5 and 2.7, even though the curves 
in Fig, 2.7 indicate that the 0.9-1. 0 sec period waves should dominate the ground 
motions at all the considered epicentral distances. 



CHAPTER HI 


SCALING RELATIONSHIPS AND AVERAGE SHAPES 
FOR NORMALIZED PSDF 


3.1 Brief Overview 

In the previous chapter, a scaling model was presented for the PSDF of 
horizontal earthquake ground motions and scaling coefficients were determined 
for three different site conditions through regression analysis. This model em- 
ployed earthquake magnitude, epicentral distance, focal depth and strong ground 
motion duration as the governing parameters. In this chapter, the scaling model 
has been presented for the PSDF normalized to fixed values of PGA and ground 
motion duration. As in the previous model, this model explicitly includes the 
effect of earthquake magnitude, epicentral distance and focal depth for three 
different geologic site conditions. In a further simplification of this model, the 
‘average’ PSDF curves corresponding to three different site conditions have been 
obtained and simple analytical forms have been proposed for these curves. 

3.2 Scaling Equation and Regression Analysis 

It is a common engineering practice to describe the severity of ground 
shaking at a site by PGA in a preliminary manner. In fact, PGA is the sole 
basis of seismic zonation in various countries around the world, even though this 
concept is slowly becoming outdated (Ttifunac (1992)). Due to this reason, the 
design spectrum curves are usually specified as conditional to a predetermined 
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value of PGA For conformity with the existing practice, it has been considered 
desirable to also obtain the scaling equation of PSDF as conditional on PGA. 
It may be recalled here that the rms value of a ground acceleration process 
may be estimated by the square root of the area under the PSDF curve for the 
acceleration process, and that by multiplying this with a suitable peak factor, 
expected PGA may be estimated. It is therefore possible to scale the ordinates 
of a given PSDF uniformly so as to correspond to an expected predetermined 
value of PGA. Further, as mentioned in the previous chapters, there are several 
definitions of strong ground motion duration available to the engineering com- 
munity, and no consensus has been reached yet on a definition that can be used 
in all the applications Due to this reason, it has been a common practice among 
researchers to consider the ground motion duration arbitarily as 20 sec in various 
PSDF-bascd studies including those on the spectrum-compatible PSDFs It has 
also been observed from the data-set (as considered in Chapter II) that most of 
the times, the estimated strong ground motion durations as per the definition 
of Trifunac and Brady {1975b) is close to 20 seC On the other hand, the appli- 
cation of the scaling model presented in Chapter II may often be inconvenient 
due to difficulty in estimating the strong ground motion duration. Keeping these 
considerations in view, an alternative scaling model is proposed to be obtained 
from the PSDF data normalized to the same PGA of 1.0^ and the same duration 
of 20 sec. It is obvious that for this model, ground motion duration ceases to be 
a governing parameter. Further, since the variation of PGA with epicentral dis- 
tance correlates very well with the attenuation characteristics in a given region, 


the attenuation function as considered in Chapter II has been omitted from this 
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alLoi native model The proposed scaling equation, can thus be expressed as 

\ogyoPSDF^(T) = 2M + bi{T)M d- b2(T)Dh 4- b^(T) + b4(T)M^ (3 1) 

where, PSDF^{T) denotes the PSDF normalized to 1.0^ PGA and 20 sec dura- 
tion, and Dh ~ V -h P?" is the hypocentral distance. It may be noted that the 
representative distance, A, as considered in the scaling relationship in Chapter 
II has been replaced by the hypocentral distance, Z)/,, This has been done for 
the simplicity in the model as the estimation of the fault-size, 5 , may not always 
be convenient in practical applications. 

For the regression analysis, the raw PSDF data-set has been prepared 
from the response spectra corresponding to the recorded accelerograms as in 
Chapter II However, here the PSDF compatible with a response spectrum has 
been computed by taking the ground motion duration equal to 20 0 sec in all 
cases. Also, the ‘average’ PSDF (as obtained by averaging the PSDFs for 
tliree response spectra) for each recorded motion has been normalized to an 
expected PGA value of 1 0^. Let this be denoted by PSDFi^i^) The re- 
gression analysis has been performed for each site condition by following the 
same method as presented in the Chapter II. Data selection has also been 
performed prior to the regression analysis in a similar way except for carry- 
ing out the bubble sorting on the basis of (IogioP5DFVr(T) — 2 M) (in place of 
(logjo PSDF{^) — 2M — 2Att{Ay M, T))) As mentioned in Chapter II, the data 
selection is aimed to minimize the ‘bias’ which might creep in the least square 
fitting due to uneven distribution of data among the magnitude classes and due 
to the relative abundance of data for some of the earthquake events. 

Let the least square estimates of the coefficients, bi(T) through 64 (P) 
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in Eq (3 1), as oVj tamed from the regression analysis at each period, T, and 
smoothed along the logioT axis, be repieseiited by Si(T) through ^^{T) respec- 
tively. Thus, the least square estimate of PSDFi^(T) may be represented by 
P^Fu{T) where 

log, „PSDF„(T) = 2M + bi{T)M + S 2 (T)Dh + SalT) + k(T)M^ . (3,2) 


This equation is considered to be applicable only in the range, Mmin ^ M < 
-Wn,ax> described in Chapter II, where 


-A^rnin — 


and 


-^max 


h{T) 

2k{T) ’ 

2 + Si(r) 
A(T) 


Thus, Eq, (3 2) may be modified to 


log,o PSDF^{T) = 

f 2M + + S2(r)Dh + S3(r) + b,{T)Ml ,, ; 

2M + bi{T)M + b2{T)Di, + k{T) + h{T)M ^ ; 

+ k(T)M„,^ + h{T)Di, + h{T) + h(T)M^ 


max ) 


(3 3) 
(3.4) 


M < M„i,n 

< M < 

^m&x ^ ^ * 


(3.5) 


Let the residual describing the deviation of the recorded normalized PSDF from 
the estimated value at period, T, be given by 


e{T) = logioPS'L>FN(r) - \Qg,^PSDFi^(T) . 


(3.6) 


Let the fluctuations in £{T) at period, T, for various PSDFs corresponding to a 
particular soil condition be assumed to be des'^ribed by a Gaussion distribution 
with mean, /i(T'), and standard deviation, cr(T). Thus, 

<T) r lyx-fi(T)y- 


<T(r)v^/-oc 2( <T(r) )] 


da; 


(3.7) 
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IS the probability of s{T) not being exceeded by the deviation at period, T. The 
actual probability, p*(e,T), that e(T') will not be exceeded has been estimated at 
each of the 91 different periods, T, from the fractions of the residues as explained 
Chapter II. The assumed probability distribution function, p(e, T), has been then 
fitted to the actual distribution function at period, T with the determination of 
the least square estimates, fi{T) and d(T), of the mean and standard deviation 
respectively for each soil condition. Each set of the estimates, fi{T) and ff(T), 
has been smoothed out along the logigT axis for estimating (normalized) PSDF 

A A A A _ 

from the smoothed estimates, l>i(T), h 2 {T), h{T) and b4(r), with a given level 
of confidence. The chi-square statistics, Kolmogorov-Smirnov 

statistics, KS(T), have also been computed at feach time period and compared 
with the respective limits for 95% level to test the ‘quality of fit’ of the normal 
distribution function described by Eq (3 7) 

3.3 Results and Discussion 

Figs. 3, 1(a) to 3 1(c) show the plots (in soild lines) of the smoothed least 
square estimates, hx{T) through b^(T)y of the regression coefficients for the allu- 
vium site, the intermediate site and the hard rock site conditions respectively. 
In each graph, dashed lines delineate the 95% confidence interval for these coef- 
ficients. Tables 3.1(a) to 3.1(c) present the (smoothed) least square estimates, 

A A 

6i(T') through 64(2^), and the (smoothed) values of the fitted normal distribution 
parameters, ii{T) and o'(T), at 12 selected periods for the three site conditions. 
The corresponding residuals as estimated by using smoothed /x(T) and d(T) val- 
ues respectively in place of /i(T') and <t{T) in Eq. (3.7), have been tabulated in 
Tables 3.2(a) to 3.2(c) for nine different confidence levels. Figs. 3 2(a) to 3,2 (c) 
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Figure 3.1(c) Regression Coefficients for Hard Rock Site 
Condition. 


central library 

* I T., KAfJPli a 






52 


Table 3 1(a) - bofust Sciuau' lislimatof! of Hogiossion GooflinfriLs and Hrisidual 
FaiauirU'i.s for Alhiviuni Silo Condition 


Peiiod, 

Least vSqua 

le Estirnatos of Rcgiession CoclTicients 

Residual Paiameteis 

T (see) 


\00ii{T) 

k{T) 

1064 ( 7 ") 

iXT) 

a(T) 

0 04 

-1 890 

-0.172 

1 800 

-0 248 

0 007 

0 304 

0 0G5 

-1.830 

-0.273 

2 170 

-0 364 

-0 014 

0 335 

0.11 

-1.800 

-0 308 

2.470 

-0 339 

-0 007 

0 340 

0 19 

-1 830 

-0 2 15 

2 440 

-0 150 

0 016 

0 324 

0.34 

-1 880 

-0 118 

2 030 

0.127 

0 031 

0.358 

0.50 

-1.900 

-0.033 

1.630 

0.276 

0 038 

0 415 

0.90 

-1.930 

0 070 

0 944 

0.456 

0.057 

0 515 

1.60 

-2.050 

0.168 

0 448 

0.675 

0 074 

0.627 

2 80 

-2.170 

0.304 

0 041 

0 853 

0 089 

0 825 

4 40 

-2.230 

0.403 

-0.361 

0 892 

0 081 

1.051 

7.50 

-2.200 

0.355 

-0.614 

0 659 

0.038 

1.271 

14.00 

-2.040 

0.036 

-0.672 

0.095 

0 006 

1 421 










Table 3.1(b) 


Lofisl S(iu<iio I^sliiriatas uf Ji 


(’girsMon Coc’fTicienLs and Ib^sidual 


Paianiolci.s for Intpi medial o Site Condition 


Pci lod , 

Regie.ss]on 

Ci)pfbcients 

Residual 

Statistics 

T (sec) 

i.(T) 

100(12(7’) 

Ut) 

l0li{T) 

,-i(T) 

a{T) 

OOd 

-2.730 

-0 277 

3.850 

0 604 

0 001 

0 267 

0 065 

-2 .360 

-0 321 

3 440 

0.179 

0 002 

0 297 

Oil 

-1 790 

-0 331 

2 400 

-0 341 

0 003 

0 302 

0 19 

-1 250 

-0 278 

1 100 

-0 717 

0 006 

0.290 

0 3d 

-1.090 

-0 123 

0 384 

-0 705 

-0 001 

0 318 

0.50 

-1 170 

0.017 

0 240 

-0.541 

-0 018 

0,386 

0 90 

-1.580 

0 200 

0 488 

-0 017 

-0 049 

0.514 

1.60 

-2,620 

0 319 

2.250 

1 080 

-0 069 

0 613 

2 80 

-3.800 

0.365 

4 670 

2 170 

-0 073 

0.708 

d.dO 

-d.740 

0 331 

6,850 

2 930 

-0.077 

0,813 

7.50 

-5.900 

0.135 

9 960 

3.700 

-0.079 

0,976 

Id.OO 

-7.680 

-0.279 

14 900 

4 960 

-0.069 

1.234 







Table 3.1(c) - Lc'fisl S(jii,u(‘ I'^limnlcs of Hcgicssion Coenicicnts and Residual 
Paiaiiictci's foi JIaid Rock Site Oondilion 


I^ciiod, 

r (see.) 

Least Squaio Estinmte.s of Rcgicssion CoefRcients 

Rc.sidual Statistics 

h ( r ) 

10062 ( 7 ’) 

63(7’) 

1064(r) 

A(T) 

a{T) 

OOd 

-1 900 

-0.173 

2 120 

-0 274 

-0.017 

0 260 

0.065 

-1 910 

-0 255 

2 490 

-0 276 

0 010 

0.338 

Oil 

-2 000 

-0 239 

2 760 

-0 051 

0.024 

0 372 

0.19 

-2 140 

-0 133 

2 600 

0 329 

0 014 

0,434 

0 34 

-2 .350 

0 207 

2.180 

0 852 

-0 014 

0 565 

0 50 

-2.460 

-0,456 

1 880 

1.110 

-0 034 

0 629 

0 90 

-2.520 

-0,836 

1.400 

1 290 

-0.063 

0.636 

1.60 

-2,570 

-0.705 

1 030 

1.340 

-0 063 

0.593 

2.80 

-2 540 

-0 339 

0 601 

1 230 

-0 036 

0.612 

4 40 

-2.520 

-0196 

0 343 

1.080 

-0 009 

0 683 

7.50 

-2.430 

-0.366 

0 241 

0.740 

0.027 

0 818 

14 00 

-2,240 

-0 745 

0 221 

0 203 

0.076 

1 046 
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1.25 


- 2.50 
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Time Period, T (sec) 

Figure 3.2(a) Residuals for p, p* = 0.1, 0.3, 0.5, 0.7 and 
0.9 in Case of Alluvium Site Condition. 
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Figure 3.2(b) Residuals for p, p* = 0.1, 0.3, 0.5, 0.7 and 
0.9 in Case of Intermediate Site Condi- 






Residues e(T) 
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Figure 3.2(c) Residuals values for p, p* = 0.1, 0.3 0.5 

0.7 and 0.9 in Case of Hard Rock Site Con- 
dition. 




Table 3-2(a) - Residual Values for Alluvium Site Condition 













Table 3.2(b) - Residual Values for Intermediate Site Condition 
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Table 3.2(c) - Residual Values for Hard Rock Site Condition 
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show the comparisons of these residuals with the actual residuals for the alluviiun 
site, intermediate site, and hard rock conditions respectively In each figure, the 
solid lines represent the actual residuals for ~ 0 1 , 0.3, 0 5, 0 7 and 0 9 wMle 
the dashed lines represent the estimated residuals for p = 0.1, 0.3, 0.5, 0.7 and 
0.9. A close agreement is observed between the actual and the estimated residU’ 
als for all the cases. To test this agreement further, let us look at the values of the 
KS{T) and x^(T) statistics as shown in Figs. 3 3(a) and (b) respectively. These 
statistics are computed at each of the 91 time periods and compared with the 
critical values of KS and tests for 95% confidence level (see the dashed lines) 
for each site condition. It is observed that the KS test passes for the alluvium 
site case at most of the periods, while in the other cases, this test passes at fewer 
periods as for the previous scaling model It is also observed from Fig 3 3(b) 
that for these site conditions, in most of the periods, the chi-square test does not 
reject the hypothesis that the error distribution is normal. For the alluvium site 
condition, however, this test fails at most of the time periods As indicated by 
these test results, the normal distribution may not be an ideal distribution How- 
ever, in view of the fact that the results for a particular site condition are passing 
one of these tests, normal distribution heis been accepted to be the reasonable 
distribution for describing the residual estimates. 

In order to illustrate how the probabilistic estimates of (normalized) 
PSDF compare with the actual PSDFs, same three example cases have been 
considered as in Chapter II. Figs. 3 4(a) to 3.4(c) show these comparisons re- 
spectively for the Morgan Hill Earthquake case, Lytle Creek Earthquake case, 
and San Fernando Earthquake case, with solid lines representing the estimated 
PSDFs corresponding to the probability levels, p = 0,1, 0 5 and 0.9 in each figure, 
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-'ieure 3.3(a) Comparison of KS Statistics at Different 
Time Periods with the 95% Cut-off Level 
for Different Site Categories. 
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Figure 3.3(b) Comparison of Chi-square Statistics at 

Different Time Periods with the 95% Cut- 
off Level for Different Site Categories. 













Figure 3.4(a) Comparison of Actual and Estimated 
PSDFs for Morgan Hill Earthquake 
Case. 



Figure 3.4(b) Comparison of Actual and Estimated 
PSDFs for Lytle Creek Earthquake 
Case. 
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Figure 3.4(c) Comparison of Actual and Estimated 
PSDFs for San Fernando Earthquake 
Case. 



and dashed lines representing the actual PSDFs. As m Chapter II, the estimated 
PSDFs have been obtained by using the least square estimates of the regression 
coefficients (along with the error estimates) at 12 time periods It may be seen 
that in all the three cases, the actual PSDFs lie witliin the 80% confidence in- 
tervals, while in the first two cases, those are in good agreement with the 50% 
confidence level estimates Here also, the 80% confidence interval is quite big, 
thus indicating the need to include more parameters in the scaling equation. It 
may however be observed that these confidence intervals are relatively narrower 
than those seen in Chapter II. Thus, it appears that the scaling relationship for 
the normalized PSDFs is slightly better. Further, it should be obvious from any 
of the Figs. 3 4(a) to 3.4(c) that the estimated PSDFs do not necessarily corre- 
spond to 1,0^ of expected PGA, even though the regression analysis was based 
on the PSDFs corresponding to this much PGA Hence, the PSDF estimated 
by using the proposed scaling relationship should be properly scaled up or down 
depending upon whether the desired PGA level is more or less than the PGA 
corresponding to the estimated PSDF. 

In order to see how the estimated PSDFs for different earthquake mag- 
nitudes compare with each other, four PSDFs curves have been estimated for M 
= 4 5, 5.5, 6.5 and 7.5 as shown in Fig. 3.5. For this, the least square estimates 
and error estimates corresponding top = 0.5 as given in Tables 3.1(a) and 3.2(a) 
respectively have been used, and other parameters have been taken as JT = 5 km 
and i? = 50 km. Fig 3.5 clearly shows the increasing dominance of long period 
waves in the PSDFs for greater magnitudes. Similar trend is observed in Fig. 3.6 
where the PSDF curves for R = 25, 50, 100 and 200 km have been compared 
with M, H and p respectively taken as 6,5, 5 km and 0 5. This is consistent with 
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i^igure 3.5 Comparison of Estimajed PSDFs for Dif- 
ferent Magnitudes in Case of Alluvium 
Site with R = 50.0 km, H = 5.0 km and 

P = 
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Figure 3.6 Comparison of Estimated PSDFs for Dif- 
ferent Epicentral Distances in Case of Al- 
luvium Site with M = 6.5, if = 5.0 km 
and p = 0,5. 





OR 


I, ho obsorvation of Tiifunac and Lee (1985b) that due to greater attenuation late 
of high frequency waves, the mo (ions recorded at large epicentral distances are 
richer in long period waves. Fig. 3 7 shows the effect of site conditions with the 
comparison of the PSDFs for alluvium, intermediate and hard rock site condi- 
tions. These PSDFs have been obtained for M = 6.5, = 50 km, H == 6 km, 

and p = 0.5. It is seen that the PSDF for alluvium site condition has the peak 
at around 0.5 sec whereas the PSDF for the hard rock condition has the jjeak 
around 0 3 sec. This is also as expected since the ground periods in the alluvium 
site conditions are generally longer compared to the site conditions with little or 
no alluvium, and hence, the shorter period motions are filtered out duidng the 
dynamic response of the alluvium as an oscillator. 

3.4 ‘Average’ PSDF Shapes 

The two models presented thus far for estimating the PSDF or normal- 
ized PSDF of earthquake ground motion are useful when the user is certain about 
the values of the governing parameters like magnitude, epicentral distance, focal 
depth etc. corresponding to the critical event In those cases, however, where 
PGA value and the geologic site characteristics are the only information available, 
non-parametric simple and crude forms of PSDF curves may be quite useful. For 
this purpose, the raw PSDF data (i e., the normalized PSDF.s corresponding to 
1.0^ expected PGA and 20 sec duration) used for the regression analysis earlier 
in this chapter has been considered and an average curve ha.s been obtained for 
each site condition. Fig. 3.8 shows the ‘average’ PSDFs for the three site con- 
ditions. These three PSDFs have also been normalized to correspond to l.Og of 
expected PGA. It may be observed in the figure that except for a sharp peak in 



Figure 3.7 Comparison of Estimated PSDPs for Dif- 
ferent Site Conditions with M = 6.5, H 
= 5.0 km and p = 0.5. 
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Figure 3.8 ‘Average’ PSDF Curves for Different Site 
Conditions. 
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the low frequency range for the alluvnim site case, the two curves corresponding 
to the alluvium and intermediate site conditions are in very good agreement. 
Both the curves show large and narrow-banded amplification of ground motions 
at high periods. The curve for hard rock site is seen to be quite different from 
these two curves as the peak in this curve has shifted towards higher frequency 
range and is much more flat, These trends are consistent with the lower val- 
ues of ground period and higher values of ground damping usually adopted in 
the Kanai-Tajimi model, as the site approaches hard rock condition. These ate 
also consistent with the trends reported in the previous regression analyses for 
FS spectra by Trifuiiac (1976b, 1979), and Trifunac and Lee (1985b) For the 
convenience of designers, the average PSDF^ curves shown in Fig. 3.8 have been 
idealized by the following simple functional forms: 

S(f) = 3700/ 0.0 </ <0,60 

= 2220 0.60 < / < 2,30 

= 2220 exp [— 0.25 (/ — 2.30)) 2.30 < / < 19,5 

= 30.12 19.5 </< 25,0 (3.8) 

for alluvium and intennediate site conditions, and 

S(f) = 428/ 0.0 < / < 2 50 

= 1070 2.50 < / < 7.0 

= 1070exp[-0.167(/ - 7.0)] 7 0<f< 25.0 (3.9) 

for hard rock conditions. In Eqs. (3 8) and (3.9), / represents the frequency in 
Hz and S{f) is the ‘average’ ground motion PSDF in cm^/sec^. It may be noted 
that both the functional forms correspond to 1 Og of expected PGA with 20 
sec ground motion duration. These forms have been compared with the actual 
■average’ curves in Figs. 3 9(a) and (b) respectively. It can be seen that the 


PSDF (m /sec ) PSDF (m /sec 



Figure 3.9(a) Comparison of Idealized and Actual 
PSDF Curves for Alluvium and Interme- 
diate Site Conditions. 



Figure 3.9(b) Comparison of Idealized and Actual 

PSDF Curves for Hard Rock Site Condi- 
tion. 
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proposed PSDF curves for the two groups of site conditions are in very good 
agreement witli the calculated ‘average’ PSDF curves. 
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CHAPTER rv 
CONCLUSIONS 


Two scaling models have been proposed m this study for the scaling of 
PSDF of horizontal ground motion m terms of earthquake magmtude, epicen- 
tral distance, focal depth These parameters can be conveniently obtained for 
practical apphcations. One of these models considers strong motion duration, £is 
defined by Trifunac and Brady (1975b), also as a govermng parameter. The other 
model has been proposed to estimate PSDF normalized to 20 sec of duration 
This model is useful where seismic hazard is specified m terms of PGA, and the 
desired PSDF has to be consistent with the specified value of PGA. Least square 
estimates of the regression coefficients and the residual estimates for different 
levels of confidence have been obtained based on the recorded motions in West- 
ern United States between 1931 and 1984 for alluvium, intermediate and hard 
rock site conditions. For obtaimng crude estimates of normalized PSDF, when 
the knowledge of PGA and geologic site conditions only is available, simple and 
convement functional relationships have also been suggested, It has been seen 
that the PSDFs estimated by proposed scaling relationships exhibit the effects 
of magnitude, epicentral distance and site conditions on the energy distribution 
as reported earlier. Also, in comparison with the duration-dependent model, the 
estimates from the model for normalized PSDF conform to a narrower range of 
variations for a given confidence interval However, both models need to consider 
greater number of source and propagation path related parameters for further 
reductions in this range. 


The jiroposed models are likely to be very useful for PSDF-based risk as- 
SGSsnient as the PSDFs estimated by these models can be assumed to coi respond 
to eqviivalcnt stationary excitations In such applications, however, transient na- 
ture of icsponse should be explicitly accounted for by using the available methods 
of analysis for SDOF and MDOF systems However, since the proposed relation- 
sliips have been obtained from the data of western United States, those should 
be applied to other seismic regions of the world only with caution. 
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